Lab 14 Homework
Isaiah Bluestein
2021-03-04
Instructions
Answer the following questions and complete the exercises in RMarkdown. Please embed all of your code and push your final work to your repository. Your final lab report should be organized, clean, and run free from errors. Remember, you must remove the # for the included code chunks to run. Be sure to add your name to the author header above. For any included plots, make sure they are clearly labeled. You are free to use any plot type that you feel best communicates the results of your analysis.
Make sure to use the formatting conventions of RMarkdown to make your report neat and clean!
Load the libraries
library(qtl)
library(qtlcharts)
library(tidyverse)
library(ggmap)1. We have a satellite collars on a number of different individuals and want to be able to quickly look at all of their recent movements at once. Please load all the data files from us_individual_collar_data and use for loop to create plots for all different individuals of the path they move on longitude and latitude.
collar <- list.files(path = "data/us_individual_collar_data", full.names = TRUE)
for (i in 1:length(data)){
data0 <- as.data.frame(read_csv(collar[i]))
print(
ggplot(data0, aes(x=long,y=lat))+
geom_path()+
geom_point()+
labs(title = collar[i], x = "longitude", y = "latitude")
)
}## Warning: Missing column names filled in: 'X1' [1]##
## ── Column specification ────────────────────────────────────────────────────────
## cols(
## X1 = col_double(),
## date = col_date(format = ""),
## collar = col_character(),
## time = col_datetime(format = ""),
## lat = col_double(),
## long = col_double()
## )
2. Please load all the data files from us_individual_collar_data and combine all data into one data frame. Create a summary to show what is the maximum and minimum of recorded data points on longitude and latitude.
all_collars <- lapply(collar, read_csv)## Warning: Missing column names filled in: 'X1' [1]##
## ── Column specification ────────────────────────────────────────────────────────
## cols(
## X1 = col_double(),
## date = col_date(format = ""),
## collar = col_character(),
## time = col_datetime(format = ""),
## lat = col_double(),
## long = col_double()
## )## Warning: Missing column names filled in: 'X1' [1]##
## ── Column specification ────────────────────────────────────────────────────────
## cols(
## X1 = col_double(),
## date = col_date(format = ""),
## collar = col_character(),
## time = col_datetime(format = ""),
## lat = col_double(),
## long = col_double()
## )## Warning: Missing column names filled in: 'X1' [1]##
## ── Column specification ────────────────────────────────────────────────────────
## cols(
## X1 = col_double(),
## date = col_date(format = ""),
## collar = col_character(),
## time = col_datetime(format = ""),
## lat = col_double(),
## long = col_double()
## )## Warning: Missing column names filled in: 'X1' [1]##
## ── Column specification ────────────────────────────────────────────────────────
## cols(
## X1 = col_double(),
## date = col_date(format = ""),
## collar = col_character(),
## time = col_datetime(format = ""),
## lat = col_double(),
## long = col_double()
## )## Warning: Missing column names filled in: 'X1' [1]##
## ── Column specification ────────────────────────────────────────────────────────
## cols(
## X1 = col_double(),
## date = col_date(format = ""),
## collar = col_character(),
## time = col_datetime(format = ""),
## lat = col_double(),
## long = col_double()
## )## Warning: Missing column names filled in: 'X1' [1]##
## ── Column specification ────────────────────────────────────────────────────────
## cols(
## X1 = col_double(),
## date = col_date(format = ""),
## collar = col_character(),
## time = col_datetime(format = ""),
## lat = col_double(),
## long = col_double()
## )## Warning: Missing column names filled in: 'X1' [1]##
## ── Column specification ────────────────────────────────────────────────────────
## cols(
## X1 = col_double(),
## date = col_date(format = ""),
## collar = col_character(),
## time = col_datetime(format = ""),
## lat = col_double(),
## long = col_double()
## )## Warning: Missing column names filled in: 'X1' [1]##
## ── Column specification ────────────────────────────────────────────────────────
## cols(
## X1 = col_double(),
## date = col_date(format = ""),
## collar = col_character(),
## time = col_datetime(format = ""),
## lat = col_double(),
## long = col_double()
## )## Warning: Missing column names filled in: 'X1' [1]##
## ── Column specification ────────────────────────────────────────────────────────
## cols(
## X1 = col_double(),
## date = col_date(format = ""),
## collar = col_character(),
## time = col_datetime(format = ""),
## lat = col_double(),
## long = col_double()
## )## Warning: Missing column names filled in: 'X1' [1]##
## ── Column specification ────────────────────────────────────────────────────────
## cols(
## X1 = col_double(),
## date = col_date(format = ""),
## collar = col_character(),
## time = col_datetime(format = ""),
## lat = col_double(),
## long = col_double()
## )full_tracking_data <- bind_rows(all_collars)
full_tracking_data## # A tibble: 240 x 6
## X1 date collar time lat long
## <dbl> <date> <chr> <dttm> <dbl> <dbl>
## 1 1 2016-02-26 A1 2016-02-26 00:00:00 37.0 -116.
## 2 2 2016-02-26 A1 2016-02-26 01:00:00 38.4 -117.
## 3 3 2016-02-26 A1 2016-02-26 02:00:00 37.9 -117.
## 4 4 2016-02-26 A1 2016-02-26 03:00:00 38.1 -115.
## 5 5 2016-02-26 A1 2016-02-26 04:00:00 35.4 -116.
## 6 6 2016-02-26 A1 2016-02-26 05:00:00 36.6 -114.
## 7 7 2016-02-26 A1 2016-02-26 06:00:00 36.9 -115.
## 8 8 2016-02-26 A1 2016-02-26 07:00:00 37.3 -115.
## 9 9 2016-02-26 A1 2016-02-26 08:00:00 35.6 -114.
## 10 10 2016-02-26 A1 2016-02-26 09:00:00 35.9 -114.
## # … with 230 more rowsfull_tracking_data %>%
summarise(max_lat = max(lat), min_lat = min(lat), max_long = max(long), min_long = min(long))## # A tibble: 1 x 4
## max_lat min_lat max_long min_long
## <dbl> <dbl> <dbl> <dbl>
## 1 41.6 26.6 -106. -123.3. Use the range of the latitude and longitude from Q2 to build an appropriate bounding box for your map and load a map from stamen in a terrain style projection and display the map. Then, build a final map that overlays the recorded path from Q1.
lat <- c(26.6116, 41.58802)
long <- c(-106.3343, -122.6017)
bbox <- make_bbox(long, lat, f = 0.5)
map <- get_map(bbox, maptype = "terrain", source = "stamen")## Map tiles by Stamen Design, under CC BY 3.0. Data by OpenStreetMap, under ODbL.ggmap(map)
for (i in 1:length(collar)){
data0 <- as.data.frame(read_csv(collar[i]))
print(
ggmap(map)+
geom_path(data = data0, aes(long, lat))+
geom_point(data = data0, aes(long, lat))+
labs(title = "Tracking Info" [i], x = "longitude", y = "latitude")
)
}## Warning: Missing column names filled in: 'X1' [1]##
## ── Column specification ────────────────────────────────────────────────────────
## cols(
## X1 = col_double(),
## date = col_date(format = ""),
## collar = col_character(),
## time = col_datetime(format = ""),
## lat = col_double(),
## long = col_double()
## )## Warning: Missing column names filled in: 'X1' [1]##
## ── Column specification ────────────────────────────────────────────────────────
## cols(
## X1 = col_double(),
## date = col_date(format = ""),
## collar = col_character(),
## time = col_datetime(format = ""),
## lat = col_double(),
## long = col_double()
## )
## Warning: Missing column names filled in: 'X1' [1]##
## ── Column specification ────────────────────────────────────────────────────────
## cols(
## X1 = col_double(),
## date = col_date(format = ""),
## collar = col_character(),
## time = col_datetime(format = ""),
## lat = col_double(),
## long = col_double()
## )
## Warning: Missing column names filled in: 'X1' [1]##
## ── Column specification ────────────────────────────────────────────────────────
## cols(
## X1 = col_double(),
## date = col_date(format = ""),
## collar = col_character(),
## time = col_datetime(format = ""),
## lat = col_double(),
## long = col_double()
## )
## Warning: Missing column names filled in: 'X1' [1]##
## ── Column specification ────────────────────────────────────────────────────────
## cols(
## X1 = col_double(),
## date = col_date(format = ""),
## collar = col_character(),
## time = col_datetime(format = ""),
## lat = col_double(),
## long = col_double()
## )
## Warning: Missing column names filled in: 'X1' [1]##
## ── Column specification ────────────────────────────────────────────────────────
## cols(
## X1 = col_double(),
## date = col_date(format = ""),
## collar = col_character(),
## time = col_datetime(format = ""),
## lat = col_double(),
## long = col_double()
## )
## Warning: Missing column names filled in: 'X1' [1]##
## ── Column specification ────────────────────────────────────────────────────────
## cols(
## X1 = col_double(),
## date = col_date(format = ""),
## collar = col_character(),
## time = col_datetime(format = ""),
## lat = col_double(),
## long = col_double()
## )
## Warning: Missing column names filled in: 'X1' [1]##
## ── Column specification ────────────────────────────────────────────────────────
## cols(
## X1 = col_double(),
## date = col_date(format = ""),
## collar = col_character(),
## time = col_datetime(format = ""),
## lat = col_double(),
## long = col_double()
## )
## Warning: Missing column names filled in: 'X1' [1]##
## ── Column specification ────────────────────────────────────────────────────────
## cols(
## X1 = col_double(),
## date = col_date(format = ""),
## collar = col_character(),
## time = col_datetime(format = ""),
## lat = col_double(),
## long = col_double()
## )
## Warning: Missing column names filled in: 'X1' [1]##
## ── Column specification ────────────────────────────────────────────────────────
## cols(
## X1 = col_double(),
## date = col_date(format = ""),
## collar = col_character(),
## time = col_datetime(format = ""),
## lat = col_double(),
## long = col_double()
## )
Load the Data
We will use the data from an experiment on hypertension in the mouse Sugiyama et al., Genomics 71:70-77, 2001
#?hyper
data(hyper)4. Create a summary of the hypertension data. How many individuals and phenotypes are included in this data set? How many gene markers and chromosomes are included in this data set? Please create a table to show the number of markers on each chromosome.
summary(hyper)## Backcross
##
## No. individuals: 250
##
## No. phenotypes: 2
## Percent phenotyped: 100 100
##
## No. chromosomes: 20
## Autosomes: 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19
## X chr: X
##
## Total markers: 174
## No. markers: 22 8 6 20 14 11 7 6 5 5 14 5 5 5 11 6 12 4 4 4
## Percent genotyped: 47.7
## Genotypes (%):
## Autosomes: BB:50.1 BA:49.9
## X chromosome: BY:53.0 AY:47.0nmar(hyper)## 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 X
## 22 8 6 20 14 11 7 6 5 5 14 5 5 5 11 6 12 4 4 45. Please make an interactive genetic map of markers for the hypertension data.
iplotMap(hyper)## Set screen size to height=700 x width=10006. Make a plot shows the pattern of missing genotype data in the hypertension dataset. Please reorder the recorded individuals according to their blood pressure phenotypes. Is there a specific pattern of missing genotype? Please explain it.
plotMissing(hyper, main = "")
plotMissing(hyper, main = "", reorder = 1)
#Individuals 50-100 are missing lots of genotype data7. Based on your answer from previous question, you probably noticed that there are gene markers without data. Please use the function drop.nullmarkers to remove markers that have no genotype data. After this, make a new summary to show the number of markers on each chromosome. How many gene markers were dropped? Where were the dropped markers located? Please use the data without nullmarkers for the following questions.
hyper_dropped <- drop.nullmarkers(hyper)
nmar(hyper_dropped)## 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 X
## 22 8 6 20 14 11 7 6 5 5 14 5 5 4 11 6 12 4 4 4nmar(hyper)## 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 X
## 22 8 6 20 14 11 7 6 5 5 14 5 5 5 11 6 12 4 4 4One marker on chomosome 14 was dropped.
8. Please conduct single-QTL analysis and create a table to give the maximum LOD score on each chromosome based on their blood pressure phenotypes. Which gene marker has the higiest LOD score? Which chromosome contains the gene marker that has the highest LOD score? Then, creates an interactive chart with LOD curves from a genome scan for all chromosomes.
hyper_dropped <- calc.genoprob(hyper_dropped, step=1)
qtl_scan <- scanone(hyper_dropped)
qtl_scan## chr pos lod
## D1Mit296 1 3.3 4.117448e-01
## c1.loc1 1 4.3 4.256930e-01
## c1.loc2 1 5.3 4.385622e-01
## c1.loc3 1 6.3 4.500490e-01
## c1.loc4 1 7.3 4.598666e-01
## c1.loc5 1 8.3 4.677428e-01
## c1.loc6 1 9.3 4.734542e-01
## c1.loc7 1 10.3 4.768268e-01
## c1.loc8 1 11.3 4.777478e-01
## c1.loc9 1 12.3 4.761722e-01
## c1.loc10 1 13.3 4.721261e-01
## c1.loc11 1 14.3 4.657072e-01
## c1.loc12 1 15.3 4.570783e-01
## c1.loc13 1 16.3 4.464742e-01
## c1.loc14 1 17.3 4.341670e-01
## c1.loc15 1 18.3 4.204707e-01
## c1.loc16 1 19.3 4.057181e-01
## D1Mit123 1 19.7 3.995971e-01
## c1.loc17 1 20.3 4.508175e-01
## c1.loc18 1 21.3 5.433096e-01
## c1.loc19 1 22.3 6.435876e-01
## c1.loc20 1 23.3 7.499236e-01
## c1.loc21 1 24.3 8.602405e-01
## c1.loc22 1 25.3 9.723147e-01
## c1.loc23 1 26.3 1.083972e+00
## c1.loc24 1 27.3 1.193255e+00
## c1.loc25 1 28.3 1.298508e+00
## c1.loc26 1 29.3 1.398416e+00
## c1.loc27 1 30.3 1.491992e+00
## c1.loc28 1 31.3 1.578553e+00
## c1.loc29 1 32.3 1.657673e+00
## D1Mit156 1 32.8 1.694372e+00
## c1.loc30 1 33.3 1.707122e+00
## c1.loc31 1 34.3 1.691829e+00
## D1Mit178 1 35.0 1.651994e+00
## c1.loc32 1 35.3 1.843815e+00
## c1.loc33 1 36.3 2.459929e+00
## D1Mit19 1 37.2 2.944148e+00
## c1.loc34 1 37.3 2.950735e+00
## c1.loc35 1 38.3 3.010586e+00
## c1.loc36 1 39.3 3.059853e+00
## c1.loc37 1 40.3 3.098968e+00
## c1.loc38 1 41.3 3.128205e+00
## D1Mit7 1 41.5 3.132878e+00
## c1.loc39 1 42.3 3.231014e+00
## c1.loc40 1 43.3 3.317505e+00
## D1Mit46 1 43.7 3.342187e+00
## D1Mit132 1 43.7 3.342187e+00
## c1.loc41 1 44.3 3.387052e+00
## c1.loc42 1 45.3 3.447292e+00
## c1.loc43 1 46.3 3.490557e+00
## c1.loc44 1 47.3 3.517783e+00
## c1.loc45 1 48.3 3.529452e+00
## D1Mit334 1 49.2 3.526704e+00
## c1.loc46 1 49.3 3.510983e+00
## c1.loc47 1 50.3 3.329680e+00
## c1.loc48 1 51.3 3.102652e+00
## c1.loc49 1 52.3 2.826315e+00
## c1.loc50 1 53.3 2.498271e+00
## c1.loc51 1 54.3 2.124402e+00
## D1Mit305 1 54.6 2.006569e+00
## c1.loc52 1 55.3 2.059839e+00
## c1.loc53 1 56.3 2.124785e+00
## c1.loc54 1 57.3 2.176021e+00
## c1.loc55 1 58.3 2.213446e+00
## c1.loc56 1 59.3 2.237311e+00
## c1.loc57 1 60.3 2.248043e+00
## c1.loc58 1 61.3 2.246155e+00
## c1.loc59 1 62.3 2.232212e+00
## c1.loc60 1 63.3 2.206845e+00
## c1.loc61 1 64.3 2.170783e+00
## D1Mit26 1 64.5 2.162361e+00
## c1.loc62 1 65.3 2.584748e+00
## c1.loc63 1 66.3 2.978416e+00
## c1.loc64 1 67.3 3.257581e+00
## D1Mit94 1 67.8 3.362918e+00
## c1.loc65 1 68.3 3.342944e+00
## c1.loc66 1 69.3 3.220749e+00
## D1Mit218 1 69.9 3.096373e+00
## c1.loc67 1 70.3 3.172136e+00
## c1.loc68 1 71.3 3.299773e+00
## c1.loc69 1 72.3 3.338664e+00
## c1.loc70 1 73.3 3.294229e+00
## c1.loc71 1 74.3 3.175861e+00
## D1Mit100 1 74.3 3.175861e+00
## c1.loc72 1 75.3 3.216347e+00
## D1Mit102 1 75.4 3.219835e+00
## c1.loc73 1 76.3 3.363839e+00
## c1.loc74 1 77.3 3.455821e+00
## c1.loc75 1 78.3 3.476983e+00
## c1.loc76 1 79.3 3.431500e+00
## c1.loc77 1 80.3 3.325261e+00
## c1.loc78 1 81.3 3.166189e+00
## D1Mit14 1 82.0 3.029099e+00
## D1Mit105 1 82.0 3.029099e+00
## D1Mit159 1 82.0 3.029217e+00
## D1Mit267 1 82.0 3.029217e+00
## c1.loc79 1 82.3 2.957400e+00
## c1.loc80 1 83.3 2.675915e+00
## c1.loc81 1 84.3 2.331015e+00
## c1.loc82 1 85.3 1.933110e+00
## c1.loc83 1 86.3 1.512487e+00
## D1Mit15 1 86.3 1.512487e+00
## c1.loc84 1 87.3 1.499329e+00
## c1.loc85 1 88.3 1.482691e+00
## c1.loc86 1 89.3 1.462261e+00
## c1.loc87 1 90.3 1.437674e+00
## c1.loc88 1 91.3 1.408508e+00
## c1.loc89 1 92.3 1.374287e+00
## c1.loc90 1 93.3 1.334515e+00
## D1Mit456 1 94.0 1.303128e+00
## c1.loc91 1 94.3 1.294538e+00
## c1.loc92 1 95.3 1.263618e+00
## c1.loc93 1 96.3 1.229213e+00
## c1.loc94 1 97.3 1.191407e+00
## c1.loc95 1 98.3 1.150334e+00
## c1.loc96 1 99.3 1.106176e+00
## c1.loc97 1 100.3 1.059167e+00
## c1.loc98 1 101.3 1.009590e+00
## c1.loc99 1 102.3 9.577744e-01
## c1.loc100 1 103.3 9.040998e-01
## c1.loc101 1 104.3 8.489877e-01
## c1.loc102 1 105.3 7.928990e-01
## c1.loc103 1 106.3 7.363254e-01
## c1.loc104 1 107.3 6.797791e-01
## c1.loc105 1 108.3 6.237795e-01
## c1.loc106 1 109.3 5.688371e-01
## c1.loc107 1 110.3 5.154374e-01
## c1.loc108 1 111.3 4.640237e-01
## c1.loc109 1 112.3 4.149806e-01
## c1.loc110 1 113.3 3.686322e-01
## c1.loc111 1 114.3 3.252182e-01
## c1.loc112 1 115.3 2.849054e-01
## D1Mit155 1 115.8 2.659435e-01
## D2Mit359 2 7.7 2.247897e-01
## c2.loc1 2 8.7 2.645578e-01
## c2.loc2 2 9.7 3.063387e-01
## c2.loc3 2 10.7 3.486135e-01
## c2.loc4 2 11.7 3.896415e-01
## c2.loc5 2 12.7 4.277313e-01
## c2.loc6 2 13.7 4.614408e-01
## c2.loc7 2 14.7 4.897297e-01
## c2.loc8 2 15.7 5.120196e-01
## c2.loc9 2 16.7 5.281822e-01
## D2Mit82 2 17.5 5.368659e-01
## c2.loc10 2 17.7 5.439595e-01
## c2.loc11 2 18.7 5.792167e-01
## c2.loc12 2 19.7 6.137385e-01
## c2.loc13 2 20.7 6.470051e-01
## c2.loc14 2 21.7 6.785253e-01
## c2.loc15 2 22.7 7.078508e-01
## c2.loc16 2 23.7 7.346046e-01
## c2.loc17 2 24.7 7.584854e-01
## c2.loc18 2 25.7 7.792732e-01
## c2.loc19 2 26.7 7.968354e-01
## c2.loc20 2 27.7 8.111152e-01
## c2.loc21 2 28.7 8.221327e-01
## D2Mit241 2 29.5 8.286552e-01
## c2.loc22 2 29.7 8.392324e-01
## c2.loc23 2 30.7 8.925159e-01
## c2.loc24 2 31.7 9.459371e-01
## c2.loc25 2 32.7 9.987163e-01
## c2.loc26 2 33.7 1.050028e+00
## c2.loc27 2 34.7 1.099035e+00
## c2.loc28 2 35.7 1.144931e+00
## c2.loc29 2 36.7 1.186969e+00
## c2.loc30 2 37.7 1.224480e+00
## c2.loc31 2 38.7 1.256902e+00
## c2.loc32 2 39.7 1.283781e+00
## c2.loc33 2 40.7 1.304779e+00
## c2.loc34 2 41.7 1.319679e+00
## c2.loc35 2 42.7 1.328380e+00
## c2.loc36 2 43.7 1.330908e+00
## c2.loc37 2 44.7 1.327412e+00
## c2.loc38 2 45.7 1.318175e+00
## c2.loc39 2 46.7 1.303611e+00
## c2.loc40 2 47.7 1.284259e+00
## c2.loc41 2 48.7 1.260772e+00
## D2Mit14 2 49.2 1.247706e+00
## c2.loc42 2 49.7 1.339822e+00
## c2.loc43 2 50.7 1.497585e+00
## c2.loc44 2 51.7 1.594565e+00
## c2.loc45 2 52.7 1.612404e+00
## c2.loc46 2 53.7 1.551728e+00
## D2Mit62 2 54.6 1.445167e+00
## c2.loc47 2 54.7 1.448669e+00
## c2.loc48 2 55.7 1.479618e+00
## c2.loc49 2 56.7 1.502134e+00
## c2.loc50 2 57.7 1.514944e+00
## c2.loc51 2 58.7 1.517074e+00
## c2.loc52 2 59.7 1.507893e+00
## c2.loc53 2 60.7 1.487131e+00
## c2.loc54 2 61.7 1.454910e+00
## c2.loc55 2 62.7 1.411760e+00
## c2.loc56 2 63.7 1.358640e+00
## c2.loc57 2 64.7 1.296937e+00
## c2.loc58 2 65.7 1.228419e+00
## c2.loc59 2 66.7 1.155125e+00
## D2Mit280 2 66.7 1.155125e+00
## c2.loc60 2 67.7 1.050643e+00
## c2.loc61 2 68.7 9.454863e-01
## c2.loc62 2 69.7 8.404592e-01
## c2.loc63 2 70.7 7.365047e-01
## c2.loc64 2 71.7 6.347035e-01
## c2.loc65 2 72.7 5.362708e-01
## c2.loc66 2 73.7 4.425293e-01
## c2.loc67 2 74.7 3.548666e-01
## c2.loc68 2 75.7 2.746698e-01
## c2.loc69 2 76.7 2.032345e-01
## c2.loc70 2 77.7 1.416981e-01
## c2.loc71 2 78.7 9.090400e-02
## c2.loc72 2 79.7 5.136724e-02
## c2.loc73 2 80.7 2.323142e-02
## c2.loc74 2 81.7 6.305447e-03
## c2.loc75 2 82.7 6.438469e-05
## c2.loc76 2 83.7 3.762446e-03
## c2.loc77 2 84.7 1.649312e-02
## D2Mit229 2 85.2 2.593693e-02
## c2.loc78 2 85.7 2.700710e-02
## c2.loc79 2 86.7 2.918845e-02
## c2.loc80 2 87.7 3.140493e-02
## c2.loc81 2 88.7 3.363052e-02
## c2.loc82 2 89.7 3.583703e-02
## c2.loc83 2 90.7 3.800070e-02
## c2.loc84 2 91.7 4.008232e-02
## c2.loc85 2 92.7 4.205817e-02
## c2.loc86 2 93.7 4.390301e-02
## c2.loc87 2 94.7 4.559538e-02
## c2.loc88 2 95.7 4.711841e-02
## c2.loc89 2 96.7 4.845748e-02
## c2.loc90 2 97.7 4.961314e-02
## D2Mit266 2 98.4 5.030940e-02
## D3Mit164 3 2.2 3.703274e-02
## c3.loc1 3 3.2 1.776899e-02
## c3.loc2 3 4.2 4.697914e-03
## c3.loc3 3 5.2 2.474810e-06
## c3.loc4 3 6.2 5.826888e-03
## c3.loc5 3 7.2 2.389910e-02
## c3.loc6 3 8.2 5.498167e-02
## c3.loc7 3 9.2 9.848260e-02
## c3.loc8 3 10.2 1.524162e-01
## c3.loc9 3 11.2 2.137576e-01
## c3.loc10 3 12.2 2.790807e-01
## c3.loc11 3 13.2 3.451588e-01
## c3.loc12 3 14.2 4.093433e-01
## c3.loc13 3 15.2 4.697214e-01
## c3.loc14 3 16.2 5.250634e-01
## c3.loc15 3 17.2 5.747443e-01
## D3Mit6 3 17.5 5.885099e-01
## c3.loc16 3 18.2 6.006449e-01
## c3.loc17 3 19.2 6.178849e-01
## c3.loc18 3 20.2 6.349102e-01
## c3.loc19 3 21.2 6.516011e-01
## c3.loc20 3 22.2 6.678342e-01
## c3.loc21 3 23.2 6.834810e-01
## c3.loc22 3 24.2 6.984131e-01
## c3.loc23 3 25.2 7.125033e-01
## c3.loc24 3 26.2 7.256272e-01
## c3.loc25 3 27.2 7.376664e-01
## c3.loc26 3 28.2 7.485101e-01
## c3.loc27 3 29.2 7.580578e-01
## c3.loc28 3 30.2 7.662212e-01
## c3.loc29 3 31.2 7.729268e-01
## c3.loc30 3 32.2 7.781173e-01
## c3.loc31 3 33.2 7.817539e-01
## c3.loc32 3 34.2 7.838177e-01
## c3.loc33 3 35.2 7.843111e-01
## c3.loc34 3 36.2 7.832582e-01
## c3.loc35 3 37.2 7.807048e-01
## D3Mit11 3 37.2 7.807048e-01
## c3.loc36 3 38.2 5.489945e-01
## c3.loc37 3 39.2 3.394395e-01
## c3.loc38 3 40.2 1.679596e-01
## c3.loc39 3 41.2 5.141111e-02
## c3.loc40 3 42.2 1.639793e-03
## c3.loc41 3 43.2 1.982390e-02
## c3.loc42 3 44.2 9.672653e-02
## D3Mit14 3 44.8 1.651760e-01
## c3.loc43 3 45.2 1.629122e-01
## c3.loc44 3 46.2 1.571982e-01
## c3.loc45 3 47.2 1.514201e-01
## c3.loc46 3 48.2 1.455899e-01
## c3.loc47 3 49.2 1.397256e-01
## c3.loc48 3 50.2 1.338445e-01
## c3.loc49 3 51.2 1.279642e-01
## c3.loc50 3 52.2 1.221028e-01
## c3.loc51 3 53.2 1.162786e-01
## c3.loc52 3 54.2 1.105095e-01
## c3.loc53 3 55.2 1.048134e-01
## c3.loc54 3 56.2 9.920511e-02
## c3.loc55 3 57.2 9.370624e-02
## D3Mit44 3 57.9 8.992874e-02
## c3.loc56 3 58.2 8.717620e-02
## c3.loc57 3 59.2 7.761290e-02
## c3.loc58 3 60.2 6.761867e-02
## c3.loc59 3 61.2 5.747075e-02
## c3.loc60 3 62.2 4.748256e-02
## c3.loc61 3 63.2 3.797641e-02
## c3.loc62 3 64.2 2.924450e-02
## c3.loc63 3 65.2 2.153903e-02
## c3.loc64 3 66.2 1.501219e-02
## D3Mit19 3 66.7 1.221712e-02
## D4Mit149 4 0.0 2.061924e+00
## c4.loc1 4 1.0 2.317346e+00
## c4.loc2 4 2.0 2.582583e+00
## c4.loc3 4 3.0 2.853421e+00
## c4.loc4 4 4.0 3.125668e+00
## c4.loc5 4 5.0 3.395519e+00
## c4.loc6 4 6.0 3.659757e+00
## c4.loc7 4 7.0 3.915797e+00
## c4.loc8 4 8.0 4.161615e+00
## c4.loc9 4 9.0 4.395645e+00
## c4.loc10 4 10.0 4.616646e+00
## c4.loc11 4 11.0 4.823589e+00
## c4.loc12 4 12.0 5.015542e+00
## c4.loc13 4 13.0 5.191582e+00
## c4.loc14 4 14.0 5.350702e+00
## D4Mit41 4 14.2 5.380397e+00
## c4.loc15 4 15.0 5.523762e+00
## c4.loc16 4 16.0 5.531748e+00
## D4Mit108 4 16.4 5.485436e+00
## c4.loc17 4 17.0 6.563385e+00
## D4Mit237 4 17.5 6.557876e+00
## c4.loc18 4 18.0 6.536449e+00
## D4Mit286 4 18.6 6.509250e+00
## c4.loc19 4 19.0 6.572621e+00
## c4.loc20 4 20.0 6.701573e+00
## c4.loc21 4 21.0 6.794140e+00
## D4Mit214 4 21.9 6.849671e+00
## c4.loc22 4 22.0 6.794676e+00
## c4.loc23 4 23.0 5.842222e+00
## D4Mit53 4 23.0 5.842222e+00
## D4Mit89 4 23.0 5.842222e+00
## c4.loc24 4 24.0 6.218547e+00
## c4.loc25 4 25.0 6.316692e+00
## D4Mit111 4 25.1 6.307900e+00
## c4.loc26 4 26.0 6.290307e+00
## c4.loc27 4 27.0 6.254650e+00
## c4.loc28 4 28.0 6.203261e+00
## D4Mit288 4 28.4 6.178436e+00
## c4.loc29 4 29.0 7.472751e+00
## D4Mit164 4 29.5 8.093661e+00
## c4.loc30 4 30.0 7.640598e+00
## D4Mit178 4 30.6 6.374922e+00
## c4.loc31 4 31.0 6.348310e+00
## D4Mit80 4 31.7 5.137589e+00
## D4Mit81 4 31.7 5.137589e+00
## c4.loc32 4 32.0 5.119870e+00
## D4Mit276 4 32.8 4.887561e+00
## c4.loc33 4 33.0 4.866401e+00
## D4Mit152 4 33.9 4.764850e+00
## c4.loc34 4 34.0 4.694565e+00
## c4.loc35 4 35.0 3.728673e+00
## D4Mit302 4 35.0 3.728673e+00
## c4.loc36 4 36.0 3.675329e+00
## c4.loc37 4 37.0 3.616906e+00
## c4.loc38 4 38.0 3.553286e+00
## c4.loc39 4 39.0 3.484368e+00
## c4.loc40 4 40.0 3.410069e+00
## c4.loc41 4 41.0 3.330329e+00
## c4.loc42 4 42.0 3.245116e+00
## c4.loc43 4 43.0 3.154439e+00
## c4.loc44 4 44.0 3.058357e+00
## c4.loc45 4 45.0 2.957005e+00
## c4.loc46 4 46.0 2.850617e+00
## c4.loc47 4 47.0 2.739565e+00
## D4Mit175 4 47.0 2.739565e+00
## c4.loc48 4 48.0 2.848913e+00
## c4.loc49 4 49.0 2.933073e+00
## c4.loc50 4 50.0 2.986540e+00
## c4.loc51 4 51.0 3.005353e+00
## c4.loc52 4 52.0 2.987290e+00
## c4.loc53 4 53.0 2.931932e+00
## c4.loc54 4 54.0 2.840757e+00
## c4.loc55 4 55.0 2.717286e+00
## c4.loc56 4 56.0 2.567186e+00
## D4Mit16 4 56.8 2.433001e+00
## c4.loc57 4 57.0 2.444104e+00
## c4.loc58 4 58.0 2.496393e+00
## c4.loc59 4 59.0 2.542522e+00
## c4.loc60 4 60.0 2.581516e+00
## c4.loc61 4 61.0 2.612546e+00
## c4.loc62 4 62.0 2.634943e+00
## c4.loc63 4 63.0 2.648206e+00
## c4.loc64 4 64.0 2.651993e+00
## c4.loc65 4 65.0 2.646116e+00
## c4.loc66 4 66.0 2.630532e+00
## c4.loc67 4 67.0 2.605338e+00
## c4.loc68 4 68.0 2.570769e+00
## c4.loc69 4 69.0 2.527198e+00
## c4.loc70 4 70.0 2.475147e+00
## c4.loc71 4 71.0 2.415286e+00
## c4.loc72 4 72.0 2.348436e+00
## c4.loc73 4 73.0 2.275560e+00
## c4.loc74 4 74.0 2.197745e+00
## D4Mit14 4 74.3 2.173612e+00
## D5Mit193 5 0.0 3.728646e-01
## c5.loc1 5 1.0 3.682948e-01
## c5.loc2 5 2.0 3.523428e-01
## c5.loc3 5 3.0 3.250162e-01
## c5.loc4 5 4.0 2.884047e-01
## c5.loc5 5 5.0 2.463322e-01
## D5Mit61 5 5.5 2.246562e-01
## c5.loc6 5 6.0 2.070219e-01
## c5.loc7 5 7.0 1.729629e-01
## c5.loc8 5 8.0 1.410057e-01
## c5.loc9 5 9.0 1.117025e-01
## c5.loc10 5 10.0 8.552311e-02
## D5Mit387 5 10.9 6.492093e-02
## c5.loc11 5 11.0 6.624613e-02
## c5.loc12 5 12.0 8.013036e-02
## c5.loc13 5 13.0 9.502738e-02
## c5.loc14 5 14.0 1.107276e-01
## D5Mit148 5 14.2 1.139451e-01
## c5.loc15 5 15.0 1.492727e-01
## D5Mit13 5 15.3 1.630031e-01
## c5.loc16 5 16.0 2.026673e-01
## c5.loc17 5 17.0 2.605245e-01
## c5.loc18 5 18.0 3.125877e-01
## D5Mit55 5 18.6 3.380438e-01
## c5.loc19 5 19.0 3.180798e-01
## c5.loc20 5 20.0 2.685648e-01
## c5.loc21 5 21.0 2.203122e-01
## c5.loc22 5 22.0 1.742967e-01
## c5.loc23 5 23.0 1.315939e-01
## c5.loc24 5 24.0 9.332448e-02
## c5.loc25 5 25.0 6.057602e-02
## c5.loc26 5 26.0 3.430530e-02
## c5.loc27 5 27.0 1.523414e-02
## c5.loc28 5 28.0 3.782072e-03
## c5.loc29 5 29.0 1.579821e-08
## c5.loc30 5 30.0 3.604356e-03
## c5.loc31 5 31.0 1.402559e-02
## c5.loc32 5 32.0 3.049083e-02
## D5Mit312 5 32.8 4.742018e-02
## c5.loc33 5 33.0 4.896144e-02
## c5.loc34 5 34.0 5.710291e-02
## c5.loc35 5 35.0 6.596335e-02
## c5.loc36 5 36.0 7.552452e-02
## c5.loc37 5 37.0 8.576019e-02
## c5.loc38 5 38.0 9.662067e-02
## c5.loc39 5 39.0 1.080533e-01
## c5.loc40 5 40.0 1.199918e-01
## c5.loc41 5 41.0 1.323605e-01
## c5.loc42 5 42.0 1.450764e-01
## c5.loc43 5 43.0 1.580515e-01
## c5.loc44 5 44.0 1.711951e-01
## c5.loc45 5 45.0 1.844164e-01
## c5.loc46 5 46.0 1.976270e-01
## c5.loc47 5 47.0 2.107422e-01
## c5.loc48 5 48.0 2.236837e-01
## c5.loc49 5 49.0 2.363800e-01
## c5.loc50 5 50.0 2.487679e-01
## c5.loc51 5 51.0 2.607913e-01
## D5Mit188 5 51.4 2.654893e-01
## c5.loc52 5 52.0 2.759096e-01
## c5.loc53 5 53.0 2.932108e-01
## c5.loc54 5 54.0 3.102936e-01
## c5.loc55 5 55.0 3.270118e-01
## c5.loc56 5 56.0 3.432283e-01
## c5.loc57 5 57.0 3.588191e-01
## c5.loc58 5 58.0 3.736763e-01
## c5.loc59 5 59.0 3.877088e-01
## c5.loc60 5 60.0 4.008459e-01
## D5Mit213 5 60.1 4.021080e-01
## c5.loc61 5 61.0 5.523188e-01
## c5.loc62 5 62.0 7.397255e-01
## c5.loc63 5 63.0 9.339516e-01
## c5.loc64 5 64.0 1.121662e+00
## c5.loc65 5 65.0 1.294591e+00
## c5.loc66 5 66.0 1.448780e+00
## D5Mit31 5 66.7 1.544770e+00
## c5.loc67 5 67.0 1.548922e+00
## c5.loc68 5 68.0 1.554072e+00
## c5.loc69 5 69.0 1.546512e+00
## c5.loc70 5 70.0 1.527286e+00
## c5.loc71 5 71.0 1.497529e+00
## c5.loc72 5 72.0 1.458415e+00
## c5.loc73 5 73.0 1.411165e+00
## D5Mit99 5 73.2 1.400852e+00
## c5.loc74 5 74.0 1.383307e+00
## D5Mit101 5 74.3 1.352388e+00
## c5.loc75 5 75.0 1.366438e+00
## c5.loc76 5 76.0 1.362578e+00
## c5.loc77 5 77.0 1.329560e+00
## D5Mit409 5 77.6 1.296112e+00
## c5.loc78 5 78.0 1.319219e+00
## c5.loc79 5 79.0 1.359789e+00
## c5.loc80 5 80.0 1.375573e+00
## c5.loc81 5 81.0 1.367060e+00
## c5.loc82 5 82.0 1.335299e+00
## D5Mit169 5 82.0 1.335299e+00
## D6Mit86 6 0.0 1.172057e-01
## c6.loc1 6 1.0 1.782575e-01
## c6.loc2 6 2.0 2.555185e-01
## c6.loc3 6 3.0 3.467135e-01
## c6.loc4 6 4.0 4.462823e-01
## c6.loc5 6 5.0 5.462297e-01
## c6.loc6 6 6.0 6.382443e-01
## c6.loc7 6 7.0 7.157566e-01
## c6.loc8 6 8.0 7.750951e-01
## c6.loc9 6 9.0 8.155407e-01
## D6Mit273 6 9.8 8.353447e-01
## c6.loc10 6 10.0 8.543155e-01
## c6.loc11 6 11.0 9.508641e-01
## c6.loc12 6 12.0 1.048968e+00
## c6.loc13 6 13.0 1.146871e+00
## c6.loc14 6 14.0 1.242921e+00
## c6.loc15 6 15.0 1.335679e+00
## c6.loc16 6 16.0 1.423980e+00
## c6.loc17 6 17.0 1.506937e+00
## c6.loc18 6 18.0 1.583916e+00
## c6.loc19 6 19.0 1.654493e+00
## c6.loc20 6 20.0 1.718414e+00
## c6.loc21 6 21.0 1.775555e+00
## D6Mit188 6 21.9 1.821166e+00
## c6.loc22 6 22.0 1.827804e+00
## c6.loc23 6 23.0 1.862040e+00
## c6.loc24 6 24.0 1.850521e+00
## c6.loc25 6 25.0 1.805372e+00
## D6Mit8 6 25.1 1.799307e+00
## c6.loc26 6 26.0 1.672506e+00
## c6.loc27 6 27.0 1.498274e+00
## c6.loc28 6 28.0 1.291016e+00
## c6.loc29 6 29.0 1.060038e+00
## D6Mit213 6 29.5 9.416038e-01
## c6.loc30 6 30.0 9.686284e-01
## c6.loc31 6 31.0 1.020105e+00
## c6.loc32 6 32.0 1.066863e+00
## c6.loc33 6 33.0 1.107428e+00
## c6.loc34 6 34.0 1.140485e+00
## c6.loc35 6 35.0 1.164955e+00
## c6.loc36 6 36.0 1.180051e+00
## c6.loc37 6 37.0 1.185327e+00
## c6.loc38 6 38.0 1.180714e+00
## c6.loc39 6 39.0 1.166563e+00
## c6.loc40 6 40.0 1.143656e+00
## D6Mit36 6 40.4 1.132296e+00
## c6.loc41 6 41.0 1.142600e+00
## c6.loc42 6 42.0 1.154658e+00
## c6.loc43 6 43.0 1.159450e+00
## c6.loc44 6 44.0 1.156153e+00
## c6.loc45 6 45.0 1.144223e+00
## c6.loc46 6 46.0 1.123453e+00
## c6.loc47 6 47.0 1.094018e+00
## c6.loc48 6 48.0 1.056490e+00
## c6.loc49 6 49.0 1.011819e+00
## c6.loc50 6 50.0 9.612791e-01
## c6.loc51 6 51.0 9.063733e-01
## D6Mit135 6 51.4 8.835489e-01
## c6.loc52 6 52.0 9.018409e-01
## c6.loc53 6 53.0 9.175342e-01
## c6.loc54 6 54.0 9.139489e-01
## c6.loc55 6 55.0 8.919032e-01
## c6.loc56 6 56.0 8.538147e-01
## D6Mit59 6 56.8 8.142006e-01
## c6.loc57 6 57.0 8.335839e-01
## c6.loc58 6 58.0 9.250587e-01
## c6.loc59 6 59.0 9.999143e-01
## c6.loc60 6 60.0 1.049059e+00
## c6.loc61 6 61.0 1.066772e+00
## c6.loc62 6 62.0 1.051934e+00
## c6.loc63 6 63.0 1.008223e+00
## D6Mit295 6 63.4 9.842685e-01
## c6.loc64 6 64.0 9.466352e-01
## c6.loc65 6 65.0 7.873967e-01
## D6Mit201 6 65.6 6.410131e-01
## c6.loc66 6 66.0 6.667164e-01
## D6Mit15 6 66.7 6.595297e-01
## D7Mit306 7 1.1 1.369747e-01
## c7.loc1 7 2.1 1.406909e-01
## c7.loc2 7 3.1 1.436216e-01
## c7.loc3 7 4.1 1.456086e-01
## c7.loc4 7 5.1 1.465171e-01
## c7.loc5 7 6.1 1.462521e-01
## c7.loc6 7 7.1 1.447699e-01
## c7.loc7 7 8.1 1.420852e-01
## c7.loc8 7 9.1 1.382735e-01
## c7.loc9 7 10.1 1.334655e-01
## c7.loc10 7 11.1 1.278367e-01
## c7.loc11 7 12.1 1.215881e-01
## c7.loc12 7 13.1 1.149397e-01
## D7Mit25 7 13.1 1.149397e-01
## c7.loc13 7 14.1 1.344947e-01
## c7.loc14 7 15.1 1.555058e-01
## c7.loc15 7 16.1 1.777400e-01
## c7.loc16 7 17.1 2.009146e-01
## c7.loc17 7 18.1 2.247041e-01
## c7.loc18 7 19.1 2.487598e-01
## c7.loc19 7 20.1 2.727286e-01
## c7.loc20 7 21.1 2.962723e-01
## c7.loc21 7 22.1 3.190823e-01
## c7.loc22 7 23.1 3.408915e-01
## c7.loc23 7 24.1 3.614809e-01
## c7.loc24 7 25.1 3.806823e-01
## c7.loc25 7 26.1 3.983802e-01
## D7Mit297 7 26.2 4.000641e-01
## c7.loc26 7 27.1 3.087214e-01
## c7.loc27 7 28.1 1.880019e-01
## D7Mit30 7 28.4 1.540588e-01
## c7.loc28 7 29.1 1.380574e-01
## c7.loc29 7 30.1 1.162435e-01
## c7.loc30 7 31.1 9.586013e-02
## c7.loc31 7 32.1 7.712150e-02
## c7.loc32 7 33.1 6.022056e-02
## c7.loc33 7 34.1 4.531666e-02
## c7.loc34 7 35.1 3.252448e-02
## c7.loc35 7 36.1 2.190708e-02
## c7.loc36 7 37.1 1.347181e-02
## D7Mit37 7 37.2 1.274700e-02
## c7.loc37 7 38.1 1.035325e-02
## c7.loc38 7 39.1 7.840182e-03
## c7.loc39 7 40.1 5.552280e-03
## c7.loc40 7 41.1 3.563860e-03
## c7.loc41 7 42.1 1.949210e-03
## c7.loc42 7 43.1 7.858521e-04
## c7.loc43 7 44.1 1.313799e-04
## c7.loc44 7 45.1 2.925572e-05
## c7.loc45 7 46.1 5.023533e-04
## c7.loc46 7 47.1 1.555770e-03
## c7.loc47 7 48.1 3.160084e-03
## c7.loc48 7 49.1 5.269858e-03
## c7.loc49 7 50.1 7.819258e-03
## c7.loc50 7 51.1 1.072929e-02
## c7.loc51 7 52.1 1.391410e-02
## c7.loc52 7 53.1 1.728721e-02
## D7Mit71 7 53.6 1.901848e-02
## c7.loc53 7 54.1 3.668836e-02
## c7.loc54 7 55.1 8.451563e-02
## D7Nds4 7 55.6 1.093732e-01
## D8Mit3 8 6.6 5.740902e-02
## c8.loc1 8 7.6 5.216647e-02
## c8.loc2 8 8.6 4.688434e-02
## c8.loc3 8 9.6 4.162552e-02
## c8.loc4 8 10.6 3.645711e-02
## c8.loc5 8 11.6 3.144780e-02
## c8.loc6 8 12.6 2.666501e-02
## c8.loc7 8 13.6 2.217189e-02
## c8.loc8 8 14.6 1.801786e-02
## c8.loc9 8 15.6 1.426458e-02
## c8.loc10 8 16.6 1.093811e-02
## c8.loc11 8 17.6 8.060957e-03
## c8.loc12 8 18.6 5.643180e-03
## c8.loc13 8 19.6 3.682956e-03
## D8Mit292 8 19.7 3.511702e-03
## c8.loc14 8 20.6 3.346407e-03
## c8.loc15 8 21.6 3.151702e-03
## c8.loc16 8 22.6 2.946934e-03
## c8.loc17 8 23.6 2.734144e-03
## c8.loc18 8 24.6 2.515750e-03
## c8.loc19 8 25.6 2.294459e-03
## c8.loc20 8 26.6 2.073154e-03
## c8.loc21 8 27.6 1.854764e-03
## c8.loc22 8 28.6 1.642133e-03
## c8.loc23 8 29.6 1.437902e-03
## c8.loc24 8 30.6 1.238807e-03
## c8.loc25 8 31.6 1.059770e-03
## c8.loc26 8 32.6 8.945253e-04
## c8.loc27 8 33.6 7.440515e-04
## D8Mit25 8 33.9 7.018864e-04
## c8.loc28 8 34.6 6.772559e-05
## c8.loc29 8 35.6 3.419713e-04
## c8.loc30 8 36.6 2.085158e-03
## c8.loc31 8 37.6 5.290088e-03
## c8.loc32 8 38.6 9.853470e-03
## c8.loc33 8 39.6 1.558406e-02
## D8Mit45 8 40.4 2.084151e-02
## c8.loc34 8 40.6 2.326602e-02
## c8.loc35 8 41.6 3.765603e-02
## c8.loc36 8 42.6 5.597442e-02
## c8.loc37 8 43.6 7.838622e-02
## c8.loc38 8 44.6 1.049773e-01
## c8.loc39 8 45.6 1.357416e-01
## c8.loc40 8 46.6 1.705860e-01
## c8.loc41 8 47.6 2.092795e-01
## c8.loc42 8 48.6 2.515240e-01
## c8.loc43 8 49.6 2.969304e-01
## c8.loc44 8 50.6 3.450475e-01
## c8.loc45 8 51.6 3.953855e-01
## c8.loc46 8 52.6 4.474384e-01
## c8.loc47 8 53.6 5.007059e-01
## c8.loc48 8 54.6 5.547094e-01
## c8.loc49 8 55.6 6.090044e-01
## c8.loc50 8 56.6 6.631885e-01
## c8.loc51 8 57.6 7.169049e-01
## c8.loc52 8 58.6 7.698438e-01
## D8Mit271 8 59.0 7.907409e-01
## c8.loc53 8 59.6 7.904592e-01
## c8.loc54 8 60.6 7.875110e-01
## c8.loc55 8 61.6 7.812066e-01
## c8.loc56 8 62.6 7.713076e-01
## c8.loc57 8 63.6 7.576405e-01
## c8.loc58 8 64.6 7.401328e-01
## c8.loc59 8 65.6 7.188156e-01
## c8.loc60 8 66.6 6.938363e-01
## c8.loc61 8 67.6 6.654643e-01
## c8.loc62 8 68.6 6.340901e-01
## c8.loc63 8 69.6 6.002165e-01
## c8.loc64 8 70.6 5.644396e-01
## c8.loc65 8 71.6 5.274153e-01
## c8.loc66 8 72.6 4.898422e-01
## c8.loc67 8 73.6 4.523829e-01
## c8.loc68 8 74.6 4.156473e-01
## D8Mit156 8 75.4 3.871433e-01
## D9Mit297 9 12.0 1.098127e-01
## c9.loc1 9 13.0 1.021163e-01
## c9.loc2 9 14.0 9.440259e-02
## c9.loc3 9 15.0 8.672402e-02
## c9.loc4 9 16.0 7.913505e-02
## c9.loc5 9 17.0 7.169101e-02
## c9.loc6 9 18.0 6.444670e-02
## c9.loc7 9 19.0 5.745495e-02
## c9.loc8 9 20.0 5.076513e-02
## c9.loc9 9 21.0 4.442178e-02
## c9.loc10 9 22.0 3.845890e-02
## c9.loc11 9 23.0 3.291850e-02
## c9.loc12 9 24.0 2.781789e-02
## D9Mit27 9 24.0 2.781789e-02
## c9.loc13 9 25.0 3.126052e-02
## c9.loc14 9 26.0 3.492258e-02
## c9.loc15 9 27.0 3.877059e-02
## c9.loc16 9 28.0 4.277742e-02
## c9.loc17 9 29.0 4.690563e-02
## c9.loc18 9 30.0 5.111212e-02
## c9.loc19 9 31.0 5.534951e-02
## c9.loc20 9 32.0 5.956779e-02
## c9.loc21 9 33.0 6.371628e-02
## c9.loc22 9 34.0 6.774561e-02
## c9.loc23 9 35.0 7.160975e-02
## c9.loc24 9 36.0 7.526771e-02
## c9.loc25 9 37.0 7.868502e-02
## c9.loc26 9 38.0 8.183473e-02
## c9.loc27 9 39.0 8.469788e-02
## c9.loc28 9 40.0 8.726358e-02
## D9Mit8 9 40.4 8.820560e-02
## c9.loc29 9 41.0 9.304230e-02
## c9.loc30 9 42.0 1.013574e-01
## c9.loc31 9 43.0 1.099021e-01
## c9.loc32 9 44.0 1.185541e-01
## c9.loc33 9 45.0 1.271831e-01
## c9.loc34 9 46.0 1.356125e-01
## c9.loc35 9 47.0 1.436838e-01
## c9.loc36 9 48.0 1.512324e-01
## c9.loc37 9 49.0 1.581052e-01
## c9.loc38 9 50.0 1.641703e-01
## c9.loc39 9 51.0 1.693250e-01
## c9.loc40 9 52.0 1.735017e-01
## c9.loc41 9 53.0 1.766709e-01
## c9.loc42 9 54.0 1.788330e-01
## c9.loc43 9 55.0 1.800494e-01
## c9.loc44 9 56.0 1.803770e-01
## D9Mit24 9 56.8 1.800599e-01
## c9.loc45 9 57.0 1.890791e-01
## c9.loc46 9 58.0 2.385554e-01
## c9.loc47 9 59.0 2.948556e-01
## c9.loc48 9 60.0 3.565589e-01
## c9.loc49 9 61.0 4.213570e-01
## c9.loc50 9 62.0 4.863050e-01
## c9.loc51 9 63.0 5.482696e-01
## c9.loc52 9 64.0 6.044208e-01
## c9.loc53 9 65.0 6.525918e-01
## c9.loc54 9 66.0 6.914535e-01
## c9.loc55 9 67.0 7.205270e-01
## c9.loc56 9 68.0 7.400833e-01
## D9Mit18 9 68.9 7.502459e-01
## D10Mit166 10 2.2 2.143123e-01
## c10.loc1 10 3.2 2.234829e-01
## c10.loc2 10 4.2 2.320672e-01
## c10.loc3 10 5.2 2.398615e-01
## c10.loc4 10 6.2 2.466793e-01
## c10.loc5 10 7.2 2.523254e-01
## c10.loc6 10 8.2 2.566576e-01
## c10.loc7 10 9.2 2.595672e-01
## c10.loc8 10 10.2 2.609931e-01
## c10.loc9 10 11.2 2.609257e-01
## c10.loc10 10 12.2 2.594076e-01
## c10.loc11 10 13.2 2.565252e-01
## c10.loc12 10 14.2 2.524221e-01
## c10.loc13 10 15.2 2.472535e-01
## D10Mit214 10 15.3 2.466849e-01
## c10.loc14 10 16.2 2.392323e-01
## c10.loc15 10 17.2 2.305273e-01
## c10.loc16 10 18.2 2.213900e-01
## c10.loc17 10 19.2 2.118219e-01
## c10.loc18 10 20.2 2.018474e-01
## c10.loc19 10 21.2 1.914949e-01
## c10.loc20 10 22.2 1.808017e-01
## c10.loc21 10 23.2 1.698145e-01
## c10.loc22 10 24.2 1.585897e-01
## c10.loc23 10 25.2 1.472143e-01
## c10.loc24 10 26.2 1.357248e-01
## c10.loc25 10 27.2 1.242196e-01
## c10.loc26 10 28.2 1.127868e-01
## c10.loc27 10 29.2 1.015190e-01
## c10.loc28 10 30.2 9.051046e-02
## c10.loc29 10 31.2 7.985533e-02
## c10.loc30 10 32.2 6.964440e-02
## c10.loc31 10 33.2 5.993827e-02
## c10.loc32 10 34.2 5.086627e-02
## c10.loc33 10 35.2 4.246558e-02
## c10.loc34 10 36.2 3.478913e-02
## c10.loc35 10 37.2 2.787630e-02
## c10.loc36 10 38.2 2.175233e-02
## c10.loc37 10 39.2 1.641164e-02
## c10.loc38 10 40.2 1.189236e-02
## c10.loc39 10 41.2 8.153760e-03
## c10.loc40 10 42.2 5.154689e-03
## c10.loc41 10 43.2 2.901991e-03
## c10.loc42 10 44.2 1.324479e-03
## c10.loc43 10 45.2 3.723625e-04
## c10.loc44 10 46.2 8.576084e-06
## c10.loc45 10 47.2 1.673225e-04
## c10.loc46 10 48.2 7.990265e-04
## c10.loc47 10 49.2 1.852895e-03
## c10.loc48 10 50.2 3.283456e-03
## D10Mit11 10 50.3 3.444413e-03
## c10.loc49 10 51.2 2.161880e-03
## c10.loc50 10 52.2 1.050129e-03
## c10.loc51 10 53.2 3.209531e-04
## c10.loc52 10 54.2 9.149230e-06
## c10.loc53 10 55.2 1.488808e-04
## c10.loc54 10 56.2 7.669277e-04
## c10.loc55 10 57.2 1.900320e-03
## c10.loc56 10 58.2 3.538624e-03
## c10.loc57 10 59.2 5.686204e-03
## c10.loc58 10 60.2 8.328484e-03
## c10.loc59 10 61.2 1.145352e-02
## c10.loc60 10 62.2 1.499818e-02
## c10.loc61 10 63.2 1.892600e-02
## c10.loc62 10 64.2 2.318246e-02
## c10.loc63 10 65.2 2.770823e-02
## c10.loc64 10 66.2 3.244171e-02
## c10.loc65 10 67.2 3.732144e-02
## c10.loc66 10 68.2 4.228826e-02
## c10.loc67 10 69.2 4.728699e-02
## D10Mit14 10 69.9 5.077499e-02
## c10.loc68 10 70.2 5.375989e-02
## c10.loc69 10 71.2 6.371340e-02
## c10.loc70 10 72.2 7.296080e-02
## c10.loc71 10 73.2 8.055670e-02
## c10.loc72 10 74.2 8.580060e-02
## c10.loc73 10 75.2 8.843148e-02
## D10Mit297 10 75.4 8.865463e-02
## D11Mit74 11 2.2 1.716721e-01
## c11.loc1 11 3.2 1.728762e-01
## c11.loc2 11 4.2 1.595362e-01
## D11Mit2 11 4.4 1.553722e-01
## D11Mit78 11 4.4 1.553722e-01
## c11.loc3 11 5.2 7.510792e-02
## c11.loc4 11 6.2 1.385575e-02
## c11.loc5 11 7.2 1.582360e-03
## c11.loc6 11 8.2 3.685685e-02
## D11Mit79 11 8.7 6.993440e-02
## D11Mit80 11 8.7 6.993440e-02
## c11.loc7 11 9.2 8.908812e-02
## c11.loc8 11 10.2 1.313214e-01
## D11Mit82 11 10.9 1.617231e-01
## c11.loc9 11 11.2 1.731660e-01
## c11.loc10 11 12.2 2.102207e-01
## D11Mit163 11 13.1 2.401204e-01
## c11.loc11 11 13.2 2.419329e-01
## c11.loc12 11 14.2 2.592164e-01
## c11.loc13 11 15.2 2.745812e-01
## c11.loc14 11 16.2 2.875719e-01
## c11.loc15 11 17.2 2.978873e-01
## D11Mit110 11 17.5 3.004353e-01
## D11Mit136 11 17.5 3.004353e-01
## c11.loc16 11 18.2 3.299695e-01
## c11.loc17 11 19.2 3.716727e-01
## D11Mit20 11 19.7 3.920578e-01
## c11.loc18 11 20.2 4.045323e-01
## c11.loc19 11 21.2 4.237471e-01
## c11.loc20 11 22.2 4.333319e-01
## c11.loc21 11 23.2 4.318544e-01
## c11.loc22 11 24.2 4.193666e-01
## D11Mit310 11 25.1 4.000344e-01
## c11.loc23 11 25.2 4.026749e-01
## c11.loc24 11 26.2 4.295613e-01
## c11.loc25 11 27.2 4.571061e-01
## c11.loc26 11 28.2 4.849386e-01
## c11.loc27 11 29.2 5.126066e-01
## c11.loc28 11 30.2 5.395947e-01
## c11.loc29 11 31.2 5.653453e-01
## c11.loc30 11 32.2 5.892724e-01
## c11.loc31 11 33.2 6.108099e-01
## c11.loc32 11 34.2 6.294357e-01
## c11.loc33 11 35.2 6.447067e-01
## c11.loc34 11 36.2 6.562852e-01
## c11.loc35 11 37.2 6.639605e-01
## c11.loc36 11 38.2 6.676603e-01
## c11.loc37 11 39.2 6.674534e-01
## c11.loc38 11 40.2 6.635416e-01
## c11.loc39 11 41.2 6.562406e-01
## c11.loc40 11 42.2 6.459620e-01
## c11.loc41 11 43.2 6.331682e-01
## D11Mit35 11 43.7 6.259821e-01
## c11.loc42 11 44.2 6.145666e-01
## c11.loc43 11 45.2 5.909498e-01
## c11.loc44 11 46.2 5.663403e-01
## c11.loc45 11 47.2 5.408207e-01
## c11.loc46 11 48.2 5.144929e-01
## c11.loc47 11 49.2 4.874785e-01
## c11.loc48 11 50.2 4.599184e-01
## c11.loc49 11 51.2 4.319724e-01
## c11.loc50 11 52.2 4.038163e-01
## c11.loc51 11 53.2 3.756395e-01
## c11.loc52 11 54.2 3.476407e-01
## c11.loc53 11 55.2 3.200226e-01
## c11.loc54 11 56.2 2.929860e-01
## c11.loc55 11 57.2 2.667236e-01
## c11.loc56 11 58.2 2.414137e-01
## c11.loc57 11 59.2 2.172142e-01
## D11Mit67 11 60.1 1.964948e-01
## c11.loc58 11 60.2 1.952491e-01
## c11.loc59 11 61.2 1.823747e-01
## c11.loc60 11 62.2 1.687886e-01
## c11.loc61 11 63.2 1.545823e-01
## c11.loc62 11 64.2 1.398819e-01
## c11.loc63 11 65.2 1.248627e-01
## c11.loc64 11 66.2 1.097067e-01
## c11.loc65 11 67.2 9.465318e-02
## c11.loc66 11 68.2 7.995953e-02
## c11.loc67 11 69.2 6.589715e-02
## c11.loc68 11 70.2 5.273640e-02
## c11.loc69 11 71.2 4.072947e-02
## c11.loc70 11 72.2 3.008148e-02
## c11.loc71 11 73.2 2.098614e-02
## c11.loc72 11 74.2 1.353153e-02
## c11.loc73 11 75.2 7.746285e-03
## c11.loc74 11 76.2 3.640601e-03
## c11.loc75 11 77.2 1.106451e-03
## c11.loc76 11 78.2 5.363517e-05
## c11.loc77 11 79.2 3.113133e-04
## c11.loc78 11 80.2 1.724442e-03
## D11Mit48 11 80.9 3.315755e-03
## D12Mit37 12 1.1 4.286562e-01
## c12.loc1 12 2.1 4.065826e-01
## c12.loc2 12 3.1 3.796842e-01
## c12.loc3 12 4.1 3.479574e-01
## c12.loc4 12 5.1 3.118241e-01
## c12.loc5 12 6.1 2.721402e-01
## c12.loc6 12 7.1 2.302699e-01
## c12.loc7 12 8.1 1.879534e-01
## c12.loc8 12 9.1 1.471220e-01
## c12.loc9 12 10.1 1.096334e-01
## c12.loc10 12 11.1 7.698911e-02
## c12.loc11 12 12.1 5.017880e-02
## c12.loc12 12 13.1 2.953213e-02
## c12.loc13 12 14.1 1.487736e-02
## c12.loc14 12 15.1 5.621659e-03
## c12.loc15 12 16.1 9.767110e-04
## D12Mit110 12 16.4 3.539764e-04
## c12.loc16 12 17.1 5.103020e-04
## c12.loc17 12 18.1 7.849582e-04
## c12.loc18 12 19.1 1.111665e-03
## c12.loc19 12 20.1 1.474904e-03
## c12.loc20 12 21.1 1.854256e-03
## c12.loc21 12 22.1 2.227583e-03
## D12Mit34 12 23.0 2.541800e-03
## c12.loc22 12 23.1 2.562578e-03
## c12.loc23 12 24.1 2.777011e-03
## c12.loc24 12 25.1 2.994878e-03
## c12.loc25 12 26.1 3.213989e-03
## c12.loc26 12 27.1 3.430221e-03
## c12.loc27 12 28.1 3.639060e-03
## c12.loc28 12 29.1 3.835849e-03
## c12.loc29 12 30.1 4.016056e-03
## c12.loc30 12 31.1 4.175568e-03
## c12.loc31 12 32.1 4.310952e-03
## c12.loc32 12 33.1 4.419689e-03
## c12.loc33 12 34.1 4.500329e-03
## c12.loc34 12 35.1 4.552563e-03
## c12.loc35 12 36.1 4.577205e-03
## c12.loc36 12 37.1 4.576063e-03
## c12.loc37 12 38.1 4.551736e-03
## c12.loc38 12 39.1 4.507350e-03
## c12.loc39 12 40.1 4.445615e-03
## D12Mit118 12 40.4 4.424810e-03
## c12.loc40 12 41.1 4.817780e-03
## c12.loc41 12 42.1 5.408506e-03
## c12.loc42 12 43.1 6.026894e-03
## c12.loc43 12 44.1 6.667743e-03
## c12.loc44 12 45.1 7.324020e-03
## c12.loc45 12 46.1 7.987881e-03
## c12.loc46 12 47.1 8.650969e-03
## c12.loc47 12 48.1 9.304755e-03
## c12.loc48 12 49.1 9.940910e-03
## c12.loc49 12 50.1 1.055167e-02
## c12.loc50 12 51.1 1.113016e-02
## c12.loc51 12 52.1 1.167065e-02
## c12.loc52 12 53.1 1.216876e-02
## c12.loc53 12 54.1 1.262148e-02
## c12.loc54 12 55.1 1.302701e-02
## c12.loc55 12 56.1 1.338554e-02
## D12Mit20 12 56.8 1.360858e-02
## D13Mit16 13 5.7 9.100472e-02
## c13.loc1 13 6.7 6.250892e-02
## c13.loc2 13 7.7 3.279096e-02
## D13Mit59 13 7.7 3.279096e-02
## c13.loc3 13 8.7 3.311359e-02
## c13.loc4 13 9.7 3.313720e-02
## c13.loc5 13 10.7 3.282075e-02
## c13.loc6 13 11.7 3.216036e-02
## c13.loc7 13 12.7 3.116943e-02
## c13.loc8 13 13.7 2.988273e-02
## c13.loc9 13 14.7 2.835235e-02
## c13.loc10 13 15.7 2.663832e-02
## D13Mit91 13 16.4 2.537065e-02
## c13.loc11 13 16.7 2.361958e-02
## c13.loc12 13 17.7 1.793911e-02
## c13.loc13 13 18.7 1.265252e-02
## c13.loc14 13 19.7 7.955132e-03
## c13.loc15 13 20.7 4.088405e-03
## c13.loc16 13 21.7 1.370349e-03
## c13.loc17 13 22.7 6.878960e-05
## c13.loc18 13 23.7 4.419798e-04
## c13.loc19 13 24.7 2.814363e-03
## c13.loc20 13 25.7 7.318752e-03
## c13.loc21 13 26.7 1.403374e-02
## c13.loc22 13 27.7 2.296820e-02
## c13.loc23 13 28.7 3.388909e-02
## c13.loc24 13 29.7 4.651540e-02
## c13.loc25 13 30.7 6.044455e-02
## c13.loc26 13 31.7 7.521457e-02
## c13.loc27 13 32.7 9.034995e-02
## c13.loc28 13 33.7 1.054049e-01
## c13.loc29 13 34.7 1.199969e-01
## c13.loc30 13 35.7 1.338268e-01
## c13.loc31 13 36.7 1.466750e-01
## c13.loc32 13 37.7 1.584439e-01
## c13.loc33 13 38.7 1.690599e-01
## c13.loc34 13 39.7 1.785267e-01
## D13Mit148 13 40.4 1.844977e-01
## c13.loc35 13 40.7 1.865074e-01
## c13.loc36 13 41.7 1.932688e-01
## c13.loc37 13 42.7 2.001162e-01
## c13.loc38 13 43.7 2.070362e-01
## c13.loc39 13 44.7 2.140152e-01
## c13.loc40 13 45.7 2.210396e-01
## c13.loc41 13 46.7 2.280960e-01
## c13.loc42 13 47.7 2.351708e-01
## c13.loc43 13 48.7 2.422509e-01
## c13.loc44 13 49.7 2.493237e-01
## c13.loc45 13 50.7 2.563771e-01
## c13.loc46 13 51.7 2.633992e-01
## c13.loc47 13 52.7 2.703792e-01
## c13.loc48 13 53.7 2.773067e-01
## c13.loc49 13 54.7 2.841721e-01
## c13.loc50 13 55.7 2.909665e-01
## c13.loc51 13 56.7 2.976816e-01
## c13.loc52 13 57.7 3.043100e-01
## c13.loc53 13 58.7 3.108450e-01
## D13Mit78 13 59.0 3.127864e-01
## D14Mit14 14 16.4 8.544850e-02
## c14.loc1 14 17.4 7.386070e-02
## c14.loc2 14 18.4 6.263623e-02
## c14.loc3 14 19.4 5.191450e-02
## c14.loc4 14 20.4 4.184067e-02
## c14.loc5 14 21.4 3.256147e-02
## c14.loc6 14 22.4 2.422012e-02
## c14.loc7 14 23.4 1.695056e-02
## c14.loc8 14 24.4 1.087150e-02
## c14.loc9 14 25.4 6.073489e-03
## c14.loc10 14 26.4 2.646345e-03
## c14.loc11 14 27.4 6.127987e-04
## c14.loc12 14 28.4 2.489631e-06
## c14.loc13 14 29.4 7.673782e-04
## c14.loc14 14 30.4 2.884577e-03
## c14.loc15 14 31.4 6.255601e-03
## c14.loc16 14 32.4 1.078880e-02
## c14.loc17 14 33.4 1.636895e-02
## c14.loc18 14 34.4 2.286995e-02
## c14.loc19 14 35.4 3.016034e-02
## D14Mit37 14 36.1 3.566251e-02
## c14.loc20 14 36.4 3.685730e-02
## c14.loc21 14 37.4 4.099915e-02
## c14.loc22 14 38.4 4.535820e-02
## c14.loc23 14 39.4 4.990819e-02
## c14.loc24 14 40.4 5.461322e-02
## c14.loc25 14 41.4 5.943103e-02
## c14.loc26 14 42.4 6.431405e-02
## c14.loc27 14 43.4 6.921073e-02
## c14.loc28 14 44.4 7.406725e-02
## c14.loc29 14 45.4 7.882937e-02
## c14.loc30 14 46.4 8.344446e-02
## c14.loc31 14 47.4 8.786344e-02
## c14.loc32 14 48.4 9.204256e-02
## c14.loc33 14 49.4 9.594486e-02
## c14.loc34 14 50.4 9.954129e-02
## c14.loc35 14 51.4 1.028113e-01
## c14.loc36 14 52.4 1.057428e-01
## D14Mit7 14 52.5 1.060171e-01
## c14.loc37 14 53.4 1.014547e-01
## c14.loc38 14 54.4 9.615812e-02
## c14.loc39 14 55.4 9.065331e-02
## c14.loc40 14 56.4 8.497877e-02
## c14.loc41 14 57.4 7.917933e-02
## c14.loc42 14 58.4 7.330508e-02
## c14.loc43 14 59.4 6.741011e-02
## c14.loc44 14 60.4 6.155081e-02
## c14.loc45 14 61.4 5.578397e-02
## c14.loc46 14 62.4 5.016482e-02
## c14.loc47 14 63.4 4.474507e-02
## c14.loc48 14 64.4 3.957111e-02
## c14.loc49 14 65.4 3.468257e-02
## c14.loc50 14 66.4 3.010736e-02
## c14.loc51 14 67.4 2.587878e-02
## D14Mit266 14 67.8 2.428618e-02
## D15Mit11 15 5.5 1.067514e+00
## D15Mit175 15 5.5 1.067514e+00
## c15.loc1 15 6.5 9.313136e-01
## c15.loc2 15 7.5 7.779486e-01
## D15Mit53 15 7.7 7.459694e-01
## c15.loc3 15 8.5 8.289017e-01
## c15.loc4 15 9.5 9.335453e-01
## c15.loc5 15 10.5 1.036974e+00
## c15.loc6 15 11.5 1.137072e+00
## c15.loc7 15 12.5 1.232173e+00
## D15Mit111 15 13.1 1.286305e+00
## c15.loc8 15 13.5 1.305523e+00
## c15.loc9 15 14.5 1.342501e+00
## c15.loc10 15 15.5 1.363015e+00
## D15Mit56 15 16.4 1.367291e+00
## D15Mit22 15 16.4 1.367291e+00
## c15.loc11 15 16.5 1.419452e+00
## c15.loc12 15 17.5 1.705369e+00
## D15Mit206 15 17.5 1.705369e+00
## D15Mit152 15 17.5 1.705369e+00
## c15.loc13 15 18.5 1.723381e+00
## c15.loc14 15 19.5 1.729778e+00
## c15.loc15 15 20.5 1.722921e+00
## c15.loc16 15 21.5 1.701518e+00
## c15.loc17 15 22.5 1.664717e+00
## c15.loc18 15 23.5 1.612210e+00
## c15.loc19 15 24.5 1.544341e+00
## c15.loc20 15 25.5 1.462256e+00
## c15.loc21 15 26.5 1.368017e+00
## c15.loc22 15 27.5 1.264627e+00
## c15.loc23 15 28.5 1.155848e+00
## c15.loc24 15 29.5 1.045753e+00
## D15Mit156 15 29.5 1.045753e+00
## c15.loc25 15 30.5 1.049079e+00
## c15.loc26 15 31.5 1.051170e+00
## c15.loc27 15 32.5 1.051899e+00
## c15.loc28 15 33.5 1.051140e+00
## c15.loc29 15 34.5 1.048775e+00
## c15.loc30 15 35.5 1.044695e+00
## c15.loc31 15 36.5 1.038801e+00
## c15.loc32 15 37.5 1.031013e+00
## c15.loc33 15 38.5 1.021267e+00
## c15.loc34 15 39.5 1.009518e+00
## c15.loc35 15 40.5 9.957464e-01
## c15.loc36 15 41.5 9.799586e-01
## c15.loc37 15 42.5 9.621878e-01
## c15.loc38 15 43.5 9.424963e-01
## c15.loc39 15 44.5 9.209757e-01
## c15.loc40 15 45.5 8.977472e-01
## c15.loc41 15 46.5 8.729597e-01
## c15.loc42 15 47.5 8.467882e-01
## c15.loc43 15 48.5 8.194303e-01
## c15.loc44 15 49.5 7.910989e-01
## c15.loc45 15 50.5 7.620288e-01
## c15.loc46 15 51.5 7.324511e-01
## c15.loc47 15 52.5 7.026007e-01
## c15.loc48 15 53.5 6.727062e-01
## c15.loc49 15 54.5 6.429827e-01
## c15.loc50 15 55.5 6.136325e-01
## D15Mit108 15 55.7 6.078228e-01
## c15.loc51 15 56.5 6.605595e-01
## c15.loc52 15 57.5 7.248587e-01
## c15.loc53 15 58.5 7.852070e-01
## c15.loc54 15 59.5 8.395240e-01
## c15.loc55 15 60.5 8.860776e-01
## c15.loc56 15 61.5 9.235321e-01
## c15.loc57 15 62.5 9.509617e-01
## D15Mit79 15 63.4 9.666624e-01
## D16Mit32 16 0.0 2.848957e-02
## c16.loc1 16 1.0 2.156787e-02
## c16.loc2 16 2.0 1.523806e-02
## c16.loc3 16 3.0 9.704954e-03
## c16.loc4 16 4.0 5.196073e-03
## c16.loc5 16 5.0 1.953900e-03
## c16.loc6 16 6.0 2.185807e-04
## c16.loc7 16 7.0 2.326060e-04
## c16.loc8 16 8.0 2.204015e-03
## c16.loc9 16 9.0 6.325185e-03
## c16.loc10 16 10.0 1.268332e-02
## c16.loc11 16 11.0 2.127774e-02
## c16.loc12 16 12.0 3.204205e-02
## c16.loc13 16 13.0 4.480478e-02
## c16.loc14 16 14.0 5.934514e-02
## c16.loc15 16 15.0 7.528067e-02
## c16.loc16 16 16.0 9.228146e-02
## c16.loc17 16 17.0 1.099784e-01
## c16.loc18 16 18.0 1.280147e-01
## c16.loc19 16 19.0 1.460477e-01
## c16.loc20 16 20.0 1.638342e-01
## c16.loc21 16 21.0 1.811162e-01
## c16.loc22 16 22.0 1.977111e-01
## c16.loc23 16 23.0 2.134854e-01
## c16.loc24 16 24.0 2.283466e-01
## c16.loc25 16 25.0 2.422518e-01
## D16Mit4 16 25.1 2.435883e-01
## c16.loc26 16 26.0 2.430528e-01
## c16.loc27 16 27.0 2.414622e-01
## c16.loc28 16 28.0 2.387976e-01
## c16.loc29 16 29.0 2.350540e-01
## c16.loc30 16 30.0 2.302517e-01
## c16.loc31 16 31.0 2.244382e-01
## D16Mit171 16 31.7 2.198069e-01
## c16.loc32 16 32.0 2.060124e-01
## D16Mit5 16 32.8 1.418847e-01
## c16.loc33 16 33.0 1.445231e-01
## c16.loc34 16 34.0 1.581482e-01
## c16.loc35 16 35.0 1.724358e-01
## c16.loc36 16 36.0 1.872940e-01
## c16.loc37 16 37.0 2.025835e-01
## c16.loc38 16 38.0 2.181574e-01
## c16.loc39 16 39.0 2.338435e-01
## c16.loc40 16 40.0 2.494546e-01
## c16.loc41 16 41.0 2.647949e-01
## c16.loc42 16 42.0 2.796683e-01
## c16.loc43 16 43.0 2.938856e-01
## c16.loc44 16 44.0 3.072724e-01
## c16.loc45 16 45.0 3.196755e-01
## c16.loc46 16 46.0 3.309683e-01
## c16.loc47 16 47.0 3.410547e-01
## c16.loc48 16 48.0 3.498711e-01
## c16.loc49 16 49.0 3.573867e-01
## c16.loc50 16 50.0 3.636025e-01
## c16.loc51 16 51.0 3.685482e-01
## D16Mit70 16 51.4 3.701823e-01
## D16Mit106 16 51.4 3.701823e-01
## D17Mit164 17 1.1 1.557465e-01
## D17Mit143 17 1.1 1.557465e-01
## D17Mit57 17 1.1 1.557464e-01
## c17.loc1 17 2.1 1.084208e-01
## D17Mit113 17 2.2 1.024779e-01
## c17.loc2 17 3.1 1.868977e-01
## D17Mit131 17 3.3 2.066626e-01
## D17Mit46 17 3.3 2.066626e-01
## c17.loc3 17 4.1 1.257509e-01
## c17.loc4 17 5.1 4.732155e-02
## D17Mit45 17 5.5 2.625086e-02
## c17.loc5 17 6.1 2.588518e-02
## D17Mit23 17 6.6 2.558274e-02
## c17.loc6 17 7.1 1.599908e-02
## c17.loc7 17 8.1 3.287372e-03
## c17.loc8 17 9.1 1.676036e-04
## c17.loc9 17 10.1 6.718869e-03
## D17Mit11 17 10.9 1.817885e-02
## c17.loc10 17 11.1 1.698206e-02
## c17.loc11 17 12.1 1.153003e-02
## c17.loc12 17 13.1 7.036932e-03
## c17.loc13 17 14.1 3.592022e-03
## c17.loc14 17 15.1 1.269086e-03
## c17.loc15 17 16.1 1.226208e-04
## c17.loc16 17 17.1 1.848966e-04
## c17.loc17 17 18.1 1.464206e-03
## c17.loc18 17 19.1 3.944487e-03
## D17Mit10 17 19.7 5.993635e-03
## c17.loc19 17 20.1 5.096804e-03
## c17.loc20 17 21.1 3.109676e-03
## c17.loc21 17 22.1 1.551131e-03
## c17.loc22 17 23.1 4.989611e-04
## c17.loc23 17 24.1 2.244376e-05
## c17.loc24 17 25.1 1.762814e-04
## c17.loc25 17 26.1 9.953380e-04
## c17.loc26 17 27.1 2.496953e-03
## c17.loc27 17 28.1 4.659855e-03
## c17.loc28 17 29.1 7.445192e-03
## c17.loc29 17 30.1 1.078867e-02
## c17.loc30 17 31.1 1.460796e-02
## c17.loc31 17 32.1 1.880918e-02
## c17.loc32 17 33.1 2.329377e-02
## D17Mit53 17 33.9 2.702012e-02
## c17.loc33 17 34.1 2.796260e-02
## c17.loc34 17 35.1 3.297911e-02
## c17.loc35 17 36.1 3.849285e-02
## c17.loc36 17 37.1 4.448119e-02
## c17.loc37 17 38.1 5.089424e-02
## c17.loc38 17 39.1 5.767034e-02
## c17.loc39 17 40.1 6.472885e-02
## c17.loc40 17 41.1 7.197504e-02
## c17.loc41 17 42.1 7.930455e-02
## c17.loc42 17 43.1 8.660875e-02
## c17.loc43 17 44.1 9.378051e-02
## c17.loc44 17 45.1 1.007196e-01
## c17.loc45 17 46.1 1.073374e-01
## c17.loc46 17 47.1 1.135600e-01
## c17.loc47 17 48.1 1.193306e-01
## c17.loc48 17 49.1 1.246100e-01
## c17.loc49 17 50.1 1.293755e-01
## D17Mit221 17 50.3 1.302640e-01
## D18Mit67 18 2.2 1.187825e-01
## c18.loc1 18 3.2 1.515378e-01
## c18.loc2 18 4.2 1.892249e-01
## c18.loc3 18 5.2 2.310017e-01
## c18.loc4 18 6.2 2.753936e-01
## c18.loc5 18 7.2 3.203454e-01
## c18.loc6 18 8.2 3.635827e-01
## c18.loc7 18 9.2 4.029606e-01
## c18.loc8 18 10.2 4.367982e-01
## c18.loc9 18 11.2 4.640661e-01
## c18.loc10 18 12.2 4.844131e-01
## c18.loc11 18 13.2 4.980918e-01
## c18.loc12 18 14.2 5.057624e-01
## D18Mit17 18 14.2 5.057624e-01
## c18.loc13 18 15.2 4.680214e-01
## c18.loc14 18 16.2 4.241576e-01
## c18.loc15 18 17.2 3.748987e-01
## c18.loc16 18 18.2 3.216382e-01
## c18.loc17 18 19.2 2.663751e-01
## c18.loc18 18 20.2 2.116295e-01
## c18.loc19 18 21.2 1.600809e-01
## c18.loc20 18 22.2 1.142001e-01
## c18.loc21 18 23.2 7.575306e-02
## c18.loc22 18 24.2 4.562701e-02
## c18.loc23 18 25.2 2.383887e-02
## c18.loc24 18 26.2 9.707361e-03
## D18Mit50 18 26.2 9.707361e-03
## c18.loc25 18 27.2 5.318071e-03
## c18.loc26 18 28.2 2.109062e-03
## c18.loc27 18 29.2 3.143456e-04
## c18.loc28 18 30.2 1.235046e-04
## c18.loc29 18 31.2 1.640058e-03
## c18.loc30 18 32.2 4.873200e-03
## c18.loc31 18 33.2 9.731596e-03
## c18.loc32 18 34.2 1.603784e-02
## c18.loc33 18 35.2 2.351771e-02
## c18.loc34 18 36.2 3.187961e-02
## c18.loc35 18 37.2 4.081412e-02
## D18Mit4 18 37.2 4.081412e-02
## D19Mit59 19 0.0 7.917868e-01
## c19.loc1 19 1.0 7.655232e-01
## c19.loc2 19 2.0 7.368230e-01
## c19.loc3 19 3.0 7.057501e-01
## c19.loc4 19 4.0 6.724239e-01
## c19.loc5 19 5.0 6.370284e-01
## c19.loc6 19 6.0 5.998006e-01
## c19.loc7 19 7.0 5.610521e-01
## c19.loc8 19 8.0 5.211546e-01
## c19.loc9 19 9.0 4.805359e-01
## c19.loc10 19 10.0 4.396617e-01
## c19.loc11 19 11.0 3.990454e-01
## c19.loc12 19 12.0 3.591838e-01
## c19.loc13 19 13.0 3.205611e-01
## c19.loc14 19 14.0 2.836172e-01
## c19.loc15 19 15.0 2.487260e-01
## c19.loc16 19 16.0 2.161772e-01
## c19.loc17 19 17.0 1.861811e-01
## D19Mit40 19 17.5 1.721773e-01
## c19.loc18 19 18.0 1.797079e-01
## c19.loc19 19 19.0 1.950936e-01
## c19.loc20 19 20.0 2.107220e-01
## c19.loc21 19 21.0 2.263123e-01
## c19.loc22 19 22.0 2.415637e-01
## c19.loc23 19 23.0 2.561581e-01
## c19.loc24 19 24.0 2.697826e-01
## c19.loc25 19 25.0 2.821495e-01
## c19.loc26 19 26.0 2.930150e-01
## c19.loc27 19 27.0 3.021933e-01
## c19.loc28 19 28.0 3.095664e-01
## c19.loc29 19 29.0 3.150869e-01
## c19.loc30 19 30.0 3.187758e-01
## c19.loc31 19 31.0 3.207154e-01
## c19.loc32 19 32.0 3.210348e-01
## D19Mit53 19 32.8 3.202411e-01
## c19.loc33 19 33.0 3.225681e-01
## c19.loc34 19 34.0 3.341613e-01
## c19.loc35 19 35.0 3.456033e-01
## c19.loc36 19 36.0 3.567559e-01
## c19.loc37 19 37.0 3.674844e-01
## c19.loc38 19 38.0 3.776462e-01
## c19.loc39 19 39.0 3.870987e-01
## c19.loc40 19 40.0 3.957143e-01
## c19.loc41 19 41.0 4.033483e-01
## c19.loc42 19 42.0 4.098943e-01
## c19.loc43 19 43.0 4.152580e-01
## c19.loc44 19 44.0 4.193671e-01
## c19.loc45 19 45.0 4.221745e-01
## c19.loc46 19 46.0 4.236595e-01
## c19.loc47 19 47.0 4.238278e-01
## c19.loc48 19 48.0 4.227113e-01
## c19.loc49 19 49.0 4.203658e-01
## c19.loc50 19 50.0 4.168624e-01
## c19.loc51 19 51.0 4.123101e-01
## c19.loc52 19 52.0 4.068088e-01
## c19.loc53 19 53.0 4.004764e-01
## c19.loc54 19 54.0 3.934364e-01
## c19.loc55 19 55.0 3.858108e-01
## D19Mit137 19 55.7 3.801940e-01
## DXMit55 X 1.1 1.556993e-03
## cX.loc1 X 2.1 5.124166e-03
## cX.loc2 X 3.1 1.096307e-02
## cX.loc3 X 4.1 1.925501e-02
## cX.loc4 X 5.1 3.012907e-02
## cX.loc5 X 6.1 4.367619e-02
## cX.loc6 X 7.1 5.992576e-02
## cX.loc7 X 8.1 7.884046e-02
## cX.loc8 X 9.1 1.003286e-01
## cX.loc9 X 10.1 1.241864e-01
## cX.loc10 X 11.1 1.501838e-01
## cX.loc11 X 12.1 1.780273e-01
## cX.loc12 X 13.1 2.073859e-01
## cX.loc13 X 14.1 2.379061e-01
## cX.loc14 X 15.1 2.692278e-01
## cX.loc15 X 16.1 3.009988e-01
## cX.loc16 X 17.1 3.328856e-01
## cX.loc17 X 18.1 3.645825e-01
## cX.loc18 X 19.1 3.958172e-01
## cX.loc19 X 20.1 4.263516e-01
## DXMit22 X 20.8 4.472058e-01
## cX.loc20 X 21.1 4.672679e-01
## cX.loc21 X 22.1 5.352306e-01
## cX.loc22 X 23.1 6.013884e-01
## cX.loc23 X 24.1 6.606192e-01
## cX.loc24 X 25.1 7.079179e-01
## cX.loc25 X 26.1 7.393047e-01
## cX.loc26 X 27.1 7.524712e-01
## cX.loc27 X 28.1 7.471670e-01
## cX.loc28 X 29.1 7.253647e-01
## DXMit16 X 29.5 7.128612e-01
## cX.loc29 X 30.1 7.419399e-01
## cX.loc30 X 31.1 7.895479e-01
## cX.loc31 X 32.1 8.349181e-01
## cX.loc32 X 33.1 8.766906e-01
## cX.loc33 X 34.1 9.136142e-01
## cX.loc34 X 35.1 9.446322e-01
## cX.loc35 X 36.1 9.689443e-01
## cX.loc36 X 37.1 9.860355e-01
## cX.loc37 X 38.1 9.956809e-01
## cX.loc38 X 39.1 9.979360e-01
## cX.loc39 X 40.1 9.931178e-01
## cX.loc40 X 41.1 9.817764e-01
## cX.loc41 X 42.1 9.646545e-01
## cX.loc42 X 43.1 9.426459e-01
## DXMit130 X 43.7 9.274964e-01Marker D4Mit164 on chromosome 4
summary(qtl_scan)## chr pos lod
## c1.loc45 1 48.3 3.529
## c2.loc45 2 52.7 1.612
## c3.loc33 3 35.2 0.784
## D4Mit164 4 29.5 8.094
## c5.loc68 5 68.0 1.554
## c6.loc23 6 23.0 1.862
## D7Mit297 7 26.2 0.400
## D8Mit271 8 59.0 0.791
## D9Mit18 9 68.9 0.750
## c10.loc8 10 10.2 0.261
## c11.loc36 11 38.2 0.668
## D12Mit37 12 1.1 0.429
## D13Mit78 13 59.0 0.313
## D14Mit7 14 52.5 0.106
## c15.loc14 15 19.5 1.730
## D16Mit70 16 51.4 0.370
## D17Mit46 17 3.3 0.207
## D18Mit17 18 14.2 0.506
## D19Mit59 19 0.0 0.792
## cX.loc38 X 39.1 0.998iplotScanone(qtl_scan)9. Based on your genome scan results, create a table which only includes those chromosomes with LOD > 1. Creates an interactive chart with LOD curves linked to estimated QTL effects for only those chromosomes with LOD > 1. Find the gene maker with the highest LOD score and describe how does the genetype of this marker influence the individual’s phenotype.
summary(qtl_scan, threshold=1)## chr pos lod
## c1.loc45 1 48.3 3.53
## c2.loc45 2 52.7 1.61
## D4Mit164 4 29.5 8.09
## c5.loc68 5 68.0 1.55
## c6.loc23 6 23.0 1.86
## c15.loc14 15 19.5 1.73plot <- iplotScanone(qtl_scan, chr = c(1, 2, 4, 5, 6, 15), hyper_dropped)
plotMarker D4Mit164 has the highest LOD, so the gene at that location will affect blood pressure phenotye the most.
10. Please save your interactive chart from Q9 as a html file hyper_iplotScanone.html and make sure your upload it to your github repository with your lab14 homework as well.
htmlwidgets::saveWidget(plot, file="hyper_iplotScanone.html")Push your final code to GitHub!
Please be sure that you check the keep md file in the knit preferences.